SOME ASYMPTOTICS FOR THE BESSEL FUNCTIONS WITH 
AN EXPLICIT ERROR TERM 
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Cn ' Abstract. We show how one can obtain an asymptotic expression for some 

special functions satisfying a second order differential equation with a very 
3 , explicit error term starting from appropriate upper bounds. We will work out 

the details for the Bessel function Jy{x) and the Airy function Ai{x) and find 

^+ ' a sharp approximation for their zeros. We also answer the question raised by 

Olcnko by showing that 
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for some explicit numerical constants ci and C2. 
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1. Introduction 
(N 

^ ' All basic formulas and asymptotic expressions for special functions we use with- 

C^ . out references can be found in [12'. To write down error terms in a compact form 

we will use 9,0i,92, ■■■■, to denote quantities with the absolute value not exceeding 
one. 

In most of the cases error terms of asymptotics of special functions are either not 
known or, at best, valid for a rather restricted rang of parameters. The following 
is a typical example of that kind (see e.g. [12, Ch. 10]). 
The Bessel function J^(a;) is defined by the series 



(1) M-) - (f ) E(-i)^-! 



{x^Ay 



5t 1 s-i^d is a solution of the following ODE: 

(2) x'^J'Jix) + xJl{x) + (a;2 - u^)J,{x) = 0. 

Theorem 1. Suppose that i/>0,x>0, w^ = ^ + ^, and let 

£1 >max(-- -,1), £2>max(---,l), 



'2 4 ' " '2 4 



then 



S^ 



"'2tH , a2 a'2tAv) 



(3) \l^ Mx)^cos{x-^,),^^^^u^—^^ 
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sin (x — uju) 
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v;^ a2i+i{v) ra Q2^2+i(^) 



, i=0 



a;^ 



2 2;2^2 + l 



where 



aiiy) 



2HI 



The assumption i^ > is not really restrictive and can be surmount by, say, 
applying the three term recurrence for J^. However if p is large or depends on x 
estimating the error term in ^ seems at least as difficult as the original task. 

In this paper we show that there is a simple way to circumvent this problem 
and to find an explicit expression for error terms which is also uniform in the 
parameters, provided one has an a priory upper bound on the absolute value of 
the considered function. In turn, in many cases such a bound may be obtained 
by using so-called Sonin's function. For Bessel and Airy functions, as well as for 
Herniite polynomials (see [2]), this can be done in a quite rutin way. For Jacobi 
and Laguerre polynomials it is a much more involved problem and the result is 
known only for oscillatory and transition regions [6], [7], [8]. It is worth noticing that 
despite the fact that it is rather a technical problem and we do have appropriate 
tools to tackle it (see e.g. Lemmas [H and S] below), one still needs a good deal 
of calculations to extend the bounds to monotonicity region. Thus, although the 
underlying idea of the method we use here is quit simple and can be applied to 
many special functions satisfying a second order differential equation, it is not 
utterly straightforward to work out the details. In this paper we will consider 
the Bessel function Jv{x) as an important example to illustrate this approach. 
We provide asymptotic expressions with an explicit error term for oscillatory and 
transition regions and also give some new estimates in the monotonicity region. As 
a corollary we derive a surprisingly accurate approximation for the Airy function 
Ai[—x)^ X > 0, and obtain an approximation with an explicit error term for its 
positive zeros, which, due to known inequalities, yields an approximation of the 
zeros of the Bessel function. 

It is worth noticing that the obtained estimates are, in a sense, best possible. 
In particular we will answer a question raised by Olenko [llj by showing that for 
v> ~h, 



(4) 



civ"^ < supa;^/^ 

x>0 



Jv{x) 



2 f TTV Tt 

COS [x 

Tra; V 2 



< C2V^ 



for some explicit numerical constants ci and C2. 

The paper is organized as follows. In the next section we describe an idea of the 
method. In section [3] we establish some upper bounds on Bessel and Airy functions 
we need in the sequel. Our main tool here will be Sonin's function. In section |4] we 
consider the error term of the standard asymptotic 



^.(^)~;i»(-T-i)' 



and prove ((4]), thus answering Olenko's question. The error term for asymptotic in 
the transition region will be derived in section [5j Finally, in section [51 using the 
approach of section[21 we establish a different sharper approximation for Bessel and 
Airy functions and their zeros. 
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2. Preliminaries 



The following approach was shortly described in [8] . We want to find an approx- 
imation of a solution of the differential equation 

(5) f" + b^{x)f{x)^0, 

in terms of some standard function F{x)^ which also satisfies a second order ODE 

Vi{F)^P2{x)F"+pi{x)F'+po{x)F^0. 
In fact, in what follows we choose F to be just cos (pix) with an appropriate function 

A possible approach to this problem is to seek for a multiplier function z{x) such 
that the differential operator 

2^2(5) = q2{x)g" + qi{x)g' + qa{x)g 

for g — g{x) — z{x)f{x) is in some sense close to I?i. For example, in the WKB-type 
approximation one chooses g{x) = \/b{x) f{x), yielding 

h' W^ ~ 9hh" 

(6) 9"~'-g' + b'gil + e{x))^0, e(x) ^ J'' . 

If e is small we can expect that g{x) is close to the solution of the equation 

h' 

" / , 1,2 n 

9o --^90 + 90 = 0, 

which is just go = McosB{x), where B{x) = J b{x)dx. 

Assume now that we have an a priori bound 15(2:) | < C. Then we can readily 
estimate the error term Ig — gol by solving ([6]) as an inhomogeneous equation. 



(7) g{x)^go-J e{t)b{t)sm{Bix)-B{t))g{t)dt, 

thus obtaining 

g = go+eC I e{t)b{t)dt. 



To derive an upper bound on |(7(a;)| we consider Sonin's function 

^'2 452 

S{x) = S{g; x)^g^ + ^.^^—^ = g- + 4^4 _ 256" + 35'^ ^"' 



then 



, _ 8b{6b'^ - 6bb'b" + b^b'") ,2 

^ y^) - /AU4 , OU2 nuuf\2 9 ■ 



(464 _,_ 3^/2 _ 255") 2 

Thus, if 46* - 266" + 36'^ > and 66'^ - 666'6" + 6^6"' > then S" > 0, and we 
obtain .g^(x) < 5 (00). 

Moreover, one can also get an upper bound on S in the following way: 

__ 6(464 + 36'^ -266") 

2(66'3 - 666'6" + 6'26"') ^ " ' 

that is 

2(66^3 _ 6j,yy/ ^ y2y//) ^ ^ ^ ^4 

■^ /'^ - (./^M , 01./?. o^,^,//^ ~ ZJZ '""" 



6(464 + 36'2 - 266") dx 464 + 36'2 - 266" 
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provided the last expression is nonnegative. Integrating, we find 



Siy).^ + m 


b"ix) 
263(2) 


1 + e{x) 


^(-)^l+l^ 


b"(y) 

263(a) 


1 + e(y) 



This shows that the envelop of g'^{x) given by S{x) is almost a constant as far as 
e(x) = o(l). 

In practically important examples the situation is somewhat more subtle as 
the coefficient b{x) may vanish. For instance for the Bessel function b{x) — 
x~^ ^x^ — i/"^ + 1/4 , and Sonin's function does not provide any information for 
the monotonicity region < a; < -^/i/^ — 1/4. Thus, one needs some supplemen- 
tary estimates to extend the bounds on \g{x)\ to this interval. Let us notice that 
although the behaviour of the solutions of (O looks less complicate in the mono- 
tonicity region, it, probably, allows only a piecewise approximation in reasonably 
simple elementary functions. 

Another rather technical problem is how to find the constants of integration in 
go- Here one either has to know, at least approximately, the value of g(x) at some 
points, e.g. at infinity, or to be able to match asymptotics in the oscillatory and 
transition regions. 

It is worth noticing that an asymptotic with an explicit error term in the transi- 
tion region, that is around a zero of b{x), can be obtained quite directly, provided 
we know a bound on |g(a;)|. Indeed, let b{a) = 0, and let d = ^ b'^{x)\ _ , then 
in a vicinity of a we can write 

b''{a + d-^/h) = <f^h + S{t), 

where 5{t) is small. The function y{t) = f{a + d^^/'^t) satisfies an Airy type ODE 

Solving it as an inhomogeneous equation one gets an explicit error term as above, 
yet facing again the problem of fixing integration constants. 

3. Upper bounds 

For the Bessel function it will be convenient to introduce the parameter /i = 
I i^^ — I; I . We also use the above notation w^ = ^ -I- ^ . 

In this section we establish some upper bounds we need in the sequel. A simplest 
inequality of this type [TH] states that for x real 

(8) JAx)< 



2Tfi/-|-l) 



For our purposes we need much more accurate estimates. To bound the Bessel 
function in the monotonicity region we will apply the following inequality given in 
[5] . We sketch a proof for self-completeness. 

Lemma 2. Let Ju{x) = x"" J,y(a;), v > — ^ , then for < a; < i^ -I- ^ , 



(9) JM. > V(2^+l)^"4x^-2z.-l > 2x 



Jt^ix) ~ 2x - 2iy + l 
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Proof. The function Jv(x) is an entire function with only real zeros satisfying the 
Laguerre inequality J'J' — JuJU > 0. Substituting here J" from the differential 
equation 

xX' + {2v + l)Jl + xj^ = 0, 
and dividing by Jy jx we obtain 

xt^[x) + {2v + \)t{x) + a; > 0, 
where t{x) = Jlj Jv Hence, 



, , , , / ^{2v+\f-'^x'^ + 2v + \ 2x 



2.T y/{2iy+iy -4.t2 + 2z/ + 1 



Since t{x) — ?► when a; — >■ 0+, whereas — ^^^- g^ > — oo, we conclude 

that 

2x 2x 

t(x) > = > 



v/(2z/+l)2-4a;2 + 2j/+l " 2zy + 1 ' 
We need the following very accurate estimate giving the value of Jy{i') 

where a = 0.09434980..., [J. 

The following theorem improves the inequality 

J^itv) <Jy{iy)fe'^^-'^'', iy>0, 0<t<l. 

given in |13j . For large v it is also stronger than the classical inequality 

/ (r) < — p-^V4(^+i) 



D 



PI p. 16], provided t > y/\n 16-2 w 0.88. 
Theorem 3. For x = tv, < i < 1 and v > 0, 

(11) JuM < Miy)^ exp ,\ / < ,,,,^,2, „^,,^ exp 



2z^ + i y 32/3r(|)t/''+i/3 ^\^2i/ + i 

Proof. By the previous lemma we have 



Jy{x) - J^ 2v + l 2v^\ 

hence 

Mty) < J.(j^)exp ; , / 
V 2i/+ 1 

which together with ([T0| yields the required result. D 

Remark 1. The function J^,j{x) = x^'^J^^x) of Lemma\^ belongs to so-called 
Pdly a- Laguerre class and satisfies the infinite series of inequalities: 

2m (2m\ 

(12) L„,(J,) = Y.{-ir+l^jU)j(^rn-,) > 0^ 
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where Li is the usual Laguerre inequality J'J' — JvJ'J > 0, (see e.g. |14j . [15j ). 
Using Lm,{J^ii) > for m > 1 leads to much more precise yet more complicated 
bounds on Jlj Jy and consequently on J^. Alternatively, one can use the inequality 
Li{J^u + ^J'l) > 0, A G M, then optimizing in A. It is worth noticing that both 
methods give an inequality .similar to (0) but in the opposite direction. Thus, one 
can use the known value o/j7iy(0) instead of Jjy{i>). 

Our main tool for bounding solutions of the second order differential equations 
will be Sonin's function. In particular, it was used by Szego to prove that for 

(13) \Mx)\ < 

V TTX 

Although he did not state this explicitly, his proof of Theorem 7.31.2. immediately 
implies 

(14) \YAx)\<\f^- 

V TTX 

His arguments go as follows: let y be a solution of the Bessel differential equation 

x'^y" + xy' + (x^ - iy'^)y = 0, 
then for |i^| < | Sonin's function 

(15) '^^^-^y' + J^.ii-^^y^ 

is increasing and inequalities (jl3|) and (jl4[) follow by calculating S{oo) from known 
asymptotics of J^ and Y,y, whereas for i/ > ^ it is decreasing for x > ^ffl, and does 
not lead, at least directly, to any explicit inequality. 

It turns out that for v > 1/2 it is more natural to deal with the function 

(16) H.(x) = |a;2-/i|i/V,(a;), 

rather than y/xJy{x). Here we will refine an inequality for the Bessel function 
obtained in [2]. First we need a bound on location of the leftmost maximum of 
H,{x). 

Lemma 4. The first positive maximum ofT-Ly{x), v > \ , is attained at a point ^ 
satisfying 



^>v^l- (2:/)-2/3 



Proof. Since obviously < ^ < /i, we can restrict ourselves to the interval (0,/i) 
and write down 

n,{x)^x''{pi-x^f'^Mx), 

where as above J,y{x) = x^" Jy{x). Then 

= n'M) = 2{f-eY/^ ^^^^ " ^')(^^^^^ + ''^ " ^'^ ^"^^^^ 

where t{x) — Jl{x) / Ji,{x) . Hence 



m 



(4j/2 - 1 _ 4^2)^ 
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and comparing this with ()11|) we obtain the inequahty 



4i/2 - 1 - 4^2 - 2 ' 

Simphfying we get 

16f-Ai2iy+l){6iy-l)^'^+{2iy~l){2v+lfi6v+l)S.^+iy{:y+l){Av'^-lf ;== p(0 > 0. 

Observe that for v > 5/3 this polynomial has the only positive zero ^07 Indeed, the 
discriminant of p(f) in ^, up to an irrelevant numerical factor, is 

v{v + l)(2zy - l)6(2zy + l)i°(108i/2 - 172zy - 5)^. 

Thus the number of positive zeros does not change for 108i^^ — 172i^ — 5 > 0, 
in particular for v > 5/3 . For v = 5/2 we obtain the following test equation 
^6 _ 21^4 ^ ;^44^2 _ 325 = 0, with the only positive zero S, w 2.14. Finally, 

p{,y) = u{'iv'^ - 2v ~ I) > 0, 

and using the substitution n = r^/2, wc find 



p{v^Jl - (2z/)-2/3 ) = p(v/r6-r4/2) = -r^r^ - lf{2r^ + 4r^ -r + 2)<0, 
hence 



e>I^Y^l-(2l/)-2/3. 

D 



Theorem 5. Let v > ^ , then for x > 0, 

(17) \x' - f,\'/VA^)\ < 

and the constant J ^ is best possible. 

Proof. For x ~ ^/Jl the result is trivial. Otherwise we shall consider three cases. 
Case 1: X > JJL. The function 

H.(a;) = (a;2-/i)i/4j,(x), 
as easy to check, satisfies the differential equation 

Ki^ - ^^—, HUx) + 4(-^ -/;);+ (6^;- A-)^ H.(x) . 0. 



Consider Sonin's function 



4x^(a;2 - /i)2 



^^") ^ ^' (") + 4(.2 :;)r+ (eS _ ^)^ ^'' (^)' 



then Ul{x) < S{x) ioi x > ^ > 0. 
One finds 



2Afix^{x^ - lj.){Ax^ +^) 



a, 



s'{x) = — ^— ^ !^ ^n'tix) > 0, 

(4(a;2 - ^)3 + (6x2 - /i)^)2 ..w - , 
hence 

\HUx)\ </\im S{x). 

Using 

w. . Ji.-i(a;) - J,y+i(a:) 
J AX) = ^ 
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and the asymptotic formula 

Jyi.^) ~ \ — cos(x ), 

V Tra; 4 

after some calculations one finds 

Since 'Hl{x) — S{x) at all local maxima the constant \/ ^ is sharp. 
Case 2: < X < y/Jl, 5 < J^ < 3. By dH]) and ly > ^ we have 



nAx) = i,-xyf^Mx)<^jy^^ 



The maximum of the last expression is attained for x = ^/v^ — v/2 , yielding 

For i^ > 5 the function T{v + 1) is increasing and it is easy to check that f{v) is 
also increasing. Hence for 2 — '^ — ^ — ^^ 0,7^ 6, we have the following estimate 

H.(x)<-^. 
r(a+ 1) 

This yields 

35/4 35/4 r^ I 

4V2r(|) 2^/2^ V TT ' 2 

2 72 [2 

H.(.)<j^<^-, 2.6<.<3. 

Case ■?.■ < x < ^, i^ > 3. Inequality (|TT|) yields 



^-^"^< 3V3r(|K+i/3 ^"PV2.+ 1 



Let r = (i^/2)^/3^ ^y Lemma [4] we can set 



X = vJl - {2v)^'^/^z = — \Jr'^ - z , r "^ < z < 1. 
This gives 

n.{x) <a{x- ^f [zr^ - 1)VV- exp (^) := A/(z), 
where A — 3^3 ,^, . We find 



/'(z) _ (l + r3+r'^-2zr'^-zV)r3 
/(z) ~ 4(l + r2)(r2-z)(zr4-l) ' 
where 

l + r'^ + r^ -2zr^ - z^r^ > (r^ - l)(l + r)(r3 - 2r'^ + r - 1) > 0, 
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for V > 19/7. Hence f{x) is increasing and 

nAx) < Am ^a(i-^J^' (r^ - l)i/V-i exp (^-^) < 

Ae-^/^(l-.-)V^exp(^)<A</f. 
This completes the proof. D 

A similar but more complicated result can be given for J^{x). 
Theorem 6. Let v > \ and x>v+ J1J3 v^l'^ , then 

x{A{x^-i.^f-ix^-lQx^v^ + ^^Y'\ ^,^ ^, ^ 2 

X V V ^ 

Proof. Let 

x%l)'^/'^(x) 

z{x) = — ^ J^ X , 

X'^ — v^ 

where 

V'(x) = 4(x^ - v^f ~ Sx** - 10x^2 + i^"*. 

First notice that '^(x) > for i/ > ^ and a; > z^ + -^i/T" i^^^'^. Indeed, the substitu- 
tions 

(19) :, = y + . + ^^.V3^ . = (2-1/3 + n)^ 

transforms V' into a polynomial in n and y with nonegative coefficients. Consider 
Sonyn's function 

Si^x) = -Ax)^^^^^^^^.'^ix). 

and its derivative 

nx) ^ ^-^(-^ - ^:)>(^)^(^) ,,.(,), 



where 



g(x) = 16(a;2 - 1.2)9 _ 4(3.2 _ ^,2)6(2^4 _ 24iy2a;2 _ 9^.4)^ 

{x^ - v^f{v^ + ^l^v^x^ + \mv^x^ + 45x*)+ 
3a;2(2iyi° + 2^v^x^ - lUiy^x^ - 66j/'*a;^ - 42iy^x^ - Sx^^), 



and 



P{x) = 

\%v^(x?- - v^f{^x^ + v^) + 8(a:2 - i.2)-''(5i.^ + 119i^V + SISj^V + \\hv^x^ + 6i/®)+ 

17!^i2 + ^\%v^^x^ + 695(^*2;'' + 492iy^x^ + 2831i^''a;^ - h^v^x^^ + ^x^'^- 

2v^{v^ + Sli^^x^ + 'iiv^x'^ + ITTi/^a;® - 18a;^). 

Applying ([T51) to P and Q we obtain polynomials with nonnegative coefficients. 

Hence, z'^{x) < S{x) for z^ > i and x > v+ 2^/3"'" '^^Z'^. Finally, using the asymptotics 

Jy{x) ~ W — cos(x -w,,), J'(a;) ^ ~\ — sin(a; -Wi,), 

V TTX V TTX 
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we obtain 

z^(x) < lini S{x) ~ — , 

a:— 7-00 TT 

and the result follows. D 

We'll need one more statement of this type for the Airy function 

M-X) = ^ (J-1/3(C) + Jl/3(C)), 



where ( — ^ — . 
Lemma 7. 

(20) {x + c)^^^Ai{-x)<—, x>-c, 

where c = 15^'^ ■ 2"^'"^. Moreover , the values of all local maxima of the function 
{x + cYi'^Ai{—x) are restricted to the interval {—i= , jg). 

Proof. Consider the function 

f{x) = {x + cy^^Ai{-x), x>-c, 
which satisfies the following differential equation: 

^''(^)-2(^^'(^)+(^+16(^)^(^)-«- 
The corresponding Sonin's function is 

S{x) = /2 + , ^ , 

^+ 16(c+a;)^ 

and 

, , , 256 cx(x + c)(x + 2c) 

b (x) ~ ^— . 

{16x{x + cy+bf 

Hence, a; = is the only maximum of S{x) for x > — c and 

fix)<f{0<^, 

where ^ = 1.12879717..., corresponds to the first maximum of f{x). Using the 
asymptotic M(— x) ^ n^^/^x^^/^, where M{x) is defined by Ai{x) — M{x) sin(/)(x), 
one finds 

lim f{x)=TT-^^^, 

a:— >oo 

and the result follows. D 

In the transition region the following inequality, which may be of independent 
interest, will be useful: 

Theorem 8. For xi,X2 > and < ly < ^, 



(21) y/x^ \J^i,{xi)J^{x2) - J^,y{x2)Jiy{xi)\ < — siuTTJ/. 

TT 
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Proof. The function 

F = F{xi,X2) = y/x^ {J-u{,Xi)Ji,{x2) - J-v{x2)Jv{xi)) 

satisfies the Bessel differential equations 

(22) A_F + 5(a;,)F = 0, i = f,2; 



where 






For a; > we consider the following majorant of F given by Sonin's function 



(^F\^ 

F^XI,X2)<S{X1,X2)^F' 



b{x2) 
By (122]) we have 



_c _ dX2 ^ ^^ ——F 

dx2 b'^{x2) \dx2 

Since 

8^2'^^'^ = ^^^' 
for 1/ < ^ , the Sonin's function is increases in 2:2. Using the asymtotics 

J,(:,) = ./X cos(x - (^!:±1)!L ) + o{x-^/'), 

\l TTX 4 

•^^ a^ = - V — sma: - ^ -^ + o x-1/2 , 

V TTX 4 

and inequalities (|13l) . ((Ti)) . on taking the limit we obtain 

F'^(x\,X2) <S{x\,X2) < lini S{x\,X2) = 

3:2—^00 

(j2(xi)-2cos7riyJ^(a;i)J_^(xi) + J2^(xi)) = 

2xi sin "Kv , ^^ , , ,^9, ,^ 4 sin -kv 
^— (J2(xi) + Y^{x^)) < ^, . 

This completes the proof. D 

To get better numerical constants the following inequalities will be useful. 
Lemma 9. For a: > 0, 

r Isintl , 2 
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Proof. We have 



sm^t , 1 1 /■°° cos2i . 1 



dt= / TTT dt < —- 



fi-)--l TT^dt. 



/o {t + x)2 2x 2 Jo (t + a;)2 2a; 

because 

f°° cos2t , ^ f^^^^ cos2t , „ 
/ 7 ^ dt = > / 7 r^ dt > 

io (t+xr f^j^ {t+x)^ 

is an alternating sum with decreasing terms. This proves 
Now we prove p4)) . Let 

I sini| 

(t + xy 
We have 

f^ xt , f°° xdt , , x + l 2 ,3, 

Next, we show that 

2 
hm f{x) = — . 

a:— foo 7r 

Using I sin(a: + 7rfc)| = sinx, < x < § , we obtain 

T(k+1) I • ,1 oo „7r • , 

dt = 






(i + a; + ttA:) 



2 



^X ''"'^„(^T7T^^'^'^^y„ V'(^)smtdt, 

where ipi^) — T'{x)/T{x) is the digamma function. Since ip'{x) = x~^ + 0{x^^) 
for a; — > oo this yields 

X f^ sint ,. „, _i X [^ smt _-^ 2 _i 



f{x) = - / dt + 0{x-') = - dt + 0{x-') = - + Oix-'). 

n Jq t + X -K Jo X TT 

Now it is enough to show that f{x) is increasing for x > j . Using the inequalities 

1 1 ,// ^ 1 1 1 



1 + X X^ 1 + X x^ {x + 1)2 



we have 



^7' (a;) = ^ (x r ^'{l^^) sintdt, 
dx \ Jq n 

V''( -)smtdt+- %}:"{ ^)sintdi = 

"■ '^ Jo '^ 

^ t + X X 1"^ t^-x 

V''( ) sinidi / ^' { )costdt> 

^ 7^ Jo I" 

2 

sin t dt 



V^ + ^ + '^ (t + x) 

' ' X TTX TTX , 

COS t dt 



t + X + TT {t + xY {t + X + Tl) 

X TTX 
h 7" TTT 1 COS tdt > 

jr/2 \'t + X + TT (t + xy 



1 
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(j + l)7r/4 

sin t dt 

JTr/i 



X + TT x^ {x + t:)^J Jq \a; + 27r (x + tt)^ 

P{x) 



COS tdt = 

Tr/2 



2Q{x)' 



where 



y'- 



TT-''P{x) = 

114688y" + 4096(259 - 2 V2 )y ^° + 1024(4216 - 77V2 )y^ + 1536(6641 - 244V2 ; 
64(242859- 16954\/2)/ +400(40011 - 5002\/2)yS + 4(2820158- 591757^2 )y^+ 
4(1332737 - U3U7V2)y^ + (1564982 - 808187V2 )y^ + 5(44796 - 40571V2 )2/2- 

3(1140 + 7109^2 )j;- 3780, 

and 

Qix) = 2/(2/ + l)2(2y+l)2(42/ + l)2(42/ + 3)2(y + 2)(2y + 3)(4y + 5)(4y + 7), y = x/n. 

It is easy to check that the polynomial P{x + | ) has only positive coefficients. 
Hence f'{x) > for x > 2/3, and f{x) < linij^^oo f{x) ~ - ■ This competes the 
proof. D 

4. Oscillatory Region 
Having at hand an upper bound on | Jiy(x)| one can estimate the difference 

/ TTX 

(25) r{x) = J — J^{x) ~ cos{x - uj^), 

in a rather elementary way. Notice that r{x) satisfies the following differential 
equation 

2-^("'-i)^''^^)' 
with the general solution of the form 

^ f^ 2^ /"°° sin(t - x) 



(26) r(a;) = cicosa; + C2sina; + W- \j-^ 1 / ttt^ Jv{t)dt. 

Now one has only to estimate the integral and to notice that as far as it is o(l), we 
have ci = C2 = by an obvious limiting argument. 
In [IT] Olenko proved the inequalities 



civ"^^^ < lim supx'^/^ 



J„{x) - \l — cos(a; - uj^] 

" TTX 



< C2V 



13/6 



and raised the question what is the best possible exponent a of j^ in this inequality. 
It turns out that that the answer a — 2 can be readily extracted from ([26|), since 
starting with a reasonably sharp approximation to the Bessel function one can 
iterate it getting more and more accurate yet more complicate approximations. 
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Theorem 10. 

(27) 
where 



Jv{x) = \ — cos(x — Wy) + Oc^x 



-3/2 



f {^" > 



a;>0, 



k' < 



2 



X> y/JI, V> 



|, 0<x<^, z/>i. 

Moreover, up to the numerical factor c the error term in {21^ is sharp. In particula 
c cannot be taken less than l/yjzn . 

Proof. To estimate the integral in ([25)1 for |;/| < ^ we apply (|24p yielding 



I^{x) 



sin(t — a;) 



J^(t)di 



< 



|sin(t-a;)| , 2 



i2 



Tree 



Thus 



2/i 



r(a;) = ci cos x + C2 sin x + i 

t:x 

For 1/ > i and a; > y/JI we use ([TT]) and the inequality arcsina; < ^ . This gives 
I sin(i — x)\ 



l^x) < 



t3/2(t2_^)l/4 



dt < 



sin (i- x)dt 
t^ 



dt 



X tv^ 



' arcsm 



Vm 



< 



y IxJJL 2x 

Hence in this case 

r(x) — ci cos X + C2 sin x + 9 
Similarly, for j/ > ^ and < a; < ,/Jl , 



2x 



I^{x) < 



x/M 



I sin(t — a;) I 

t3/2(i2_^)l/4 



dt 



\/m 



I sin(t — a:) I 

, t3/2(/i-t2)l/' 



■dt < 



1 /¥ 
2Vm 



v^ |sin(t-a;)| 1 /¥ 

., i3/2(^_i2)l/4«^- 2 Y/i' 



t2 



x/m 



dt 



Wli - 1^ 



1 /F 
2Vm 



2x 



2^ 



/iX X 'ZX Xy/JI 

An elementary investigation shows that the maximum of the function 

x(VM-.)ln^ 

X 

is attained for x = d^/JI , d — 0.314711..., and does not exceed 2///5. Hence, 

I.(x) < ^^ < 
and 



r{x) — C\ cos X + C2 sin x + 



7x ' 



7x 
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Finally, since for a; — >■ oo, 



J^{x) = \ — cos{x ~ Lu^) + 0{x '^/^), 

V TTX 

that is lima;^oo r(x) — 0, we conclude that ci — C2 — and (|27|) follows. 

Let us show now that up to the numerical factor c the error term in (j27p is sharp. 
By (HH) and ^ we have 



n^{x) =x^/^ 



Jv{x) - \ — cos(a; - u^) 
nx 



lix 



sm(t — x) 



t3/2 



J^{t)dt 



— iix 



sin(i — x) cos(i — uj^) 



t^ 



dt 



-6cfix 



sin(t — x) 



P 



dt 



:= \l - fj.x\Ii\+l2- 



Here 



I2 — Ocfix 



dt Ocfj, 



To bound /i we introduce two auxiliary functions / and g (see [HI Ch. 6]), defined 

by 

TT 

Si{x) = — — f{x) COS a; — g{x) sin a:, 



Ci{x) — f{x) sin a; — ^(a;) cos a;, 



with the asymptotics 



/(a;) = - + 0(a; ^), g(a;) = ^ + 0(a; ''), x -^ 00. 



Calculations yield 

/i = 5*1(22) sin(x + cji,) + Ci{2z) cos(a; + uj„) - 

f{2z) sin(2i + a: — cui,) — g{2z) cos(2t + a; — uj^) 
sin{LL>i, — x) 



sintcos(2: — Wj/) 



sinicos(z — lUi,) 



t=o 



2x 



+ 0{x ^), z^x + t. 



Hence 



M 



K^{x) — — \sin{LJiy — x)\ + 0{x ^), 



/27r 



and the result follows. D 

Applying (j27p to the Airy function yields the following approximation: 
Corollary 1. For x > 0, 



(28) 



Ai{-x) 



cos(C - f ) 



/^a;i/4 6\/3 7r3/2a;7/4 ' 
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5. Transition Region 

Our estimates in the transition region are based on the following simple obser- 
vation. 

Lemma 11. Let for a real constant q, 

y"{x) + q xy{x) — u{x). 

Then for a; > 0, 

(29) y(x) = V^ j^ci J_i/3(^i^— ) + C2 Ji/3(^^) j + 

9q-^x-^/^ r \u{t)\t-^'^dt, 
Jo 
provided the integral exists. 

Proof. Let yi and 2/2 be two linearly independent solutions of the homogeneous 
equation 

y"{x)+q^xy{x)=Q, 

and let 

rj, ,N ^ Vi{i)V2{x) -~yi{x)y2{t) 

^'' yiiMit) - y'iit)y2{t) - 

Then 

y{x) = ciyi(x) + 022/2(2;) + / U{x, t)u{t)dt, 
and choosing 

2qx^ 2grV2 

2/1 (x) = y/xJ_i/3( ), y2(x) = VxJl/3[ ), 

we find 

yi(%2W-y'i(%2(t) = ^. 

Applying (PT|) we obtain 



... ^ . /^^. .2qt3/2 9g:j3/2 ,2ga;3/\ ,2gt3/2 \ 

f/(x,t) - 2n^- (^J_i/3( JL_)Ji/3(Ji-_) - J_^/3(Jl__)J^/3(JL_)j 

Hence 



U{x,t)u{t)dt 




<q-\-^'^ / \u{t)\t-^'Ut, 
Jo 



and the result follows. D 

To find the constant of integration we use the value of J^{x) at two points: x — i^ 
and at the first zero x = ji,i. The following bounds are known [16] (see [3] for a 
review of recent results), 

(30) j-^^ = ^ + 2-1/3^,^1/3 + 02 3:2_2^^_,/3^ ^^^^ 

where Ug is s*'' positive zero of the Airy function Ai{—x). 
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In the sequel we set 7 = 2^^l'^a\ — 1.855757... Thus, we can write 

(31) U = v + iv^i^^e^^-^v-^i^. 

Lemma 12. 

(32) j^(^ + V/3)< 7^ ^>1 

\Sv 2 

Proof. Let j — j„i, by |3T|) we have 

J,(. + 7.V3) ^ _02 ^ ,-1/3 J^(, + ^,1/3 + ^2 ^ ,-1/3), 

Settmg 



by (fTS]) and 7 > -^i/F" we obtain 






where 

</)(e, S) = ^^ (4(8.5^ - 3) + 16{3S^ - 2)e + 4(6(53 _ 7)e^ + A{d'^ - 'S)e^ - 3e^)"'^^ 

Notice that (f){e, S) < 0(e, 7) since ^ > 7 and 

^0(.,^) = ^^0^(e,5)<O. 



Moreover, 



6„-y^ = 71.-2/3 + ^ ,-4/3 < 6^ 



whereas 

^</)(e,7)<0, 0<e<7. 
Thus, (/)(e,5) < (f>{0,j), and we get 

and the result follows. D 

Applying Lemma [11] to Ju we obtain 
Theorem 13. Let v > ^ , and let x = v + v^/^z, then 

,33) ..M = ^|^A.(-2./3.)..f°;"";f''. 

Vj/2/3 + 2 2i/2/3Vj/2/3 + 2 

Proof. Consider the function 

y{z) = \/v + i/i/3z J^iv + i^^/^z), 
which, as easy to check, satisfies the differential equation 

„^ ^ o . ^ 8z3 + 12iy2/322 _ 1 
y (z) + 2zy(z) = 4(^2/3 + ^)2 ^ W' 
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Lemmas [IT] and [5] yield 

(34) j/(z) = Vi(ciJ_i/3(V2C)+C2Ji/3(V2C)) +0i?(z), 

where C = ^ — , and 



z'^'R{z)= ^ 



4.V2J0 



8t^ + 12j/2/3t2 - 1 

(i.2/3 + t)2tl/4 y^^> 



dt < 



T7^-(^ -^) ' y — dt 



4V2 Jo (y'^'^ + i)3/2ii/4 

Z/V6 /- |8i3 + 12^.2/3^2 „ 1 



2v^ Jo iy'^''^ + 03/2(41.2/3^2 4^ 8^4/3i + 1)1/4^1/ 



dt < 



^ r r (3t/2/3 + 2t)t3/2dt /•" d±_ 

f^'^'^^Wo {iy^/^ + t)^/H2l,V3+ty/i^J^ 4Vt (1^2/3 +t)3/2 (21.2/3 +i) 1/4 /"< 

1 r r (3i/2/3 + 2t)i3/2dt [^ dt 

max ■ 



V^ [Jo (i.2/3+f)3/2(2j.2/3 + i)l/4'_yg 4-21/4^,7/6^ 

For i > the function 

3t/2/3 + 2t 

(1/2/3 +i)3/2(2i.2/3_^i)l/4 

is decreasing in t and therefore is less than 3 • 2^^^'^v^^/'^ , its value at i = 0. Thus, 
we obtain 

1 ( f^ 3i'^/2 /z 

(35) i?(z) < ^=-^ max jj^ ^^^^ dt, 



1/4. a; ' 25/4j,7/6 



max<^ - z^^'^,1^ ^/3 



25/4 0Fi; ]^ 5 

It is left to find the constants ci,C2. For x — u, that is for z = 0, we have 
R{z) — 0, and comparing ([M)) with PH)) we get 

2/(0) = Inn Vl(ciJ_i/3(V2C) +C2Ji/3(y2C)) = ^i 1™, V^ ' ^1/6^2 )^ = 
3i/3ci ^^,, 21/307 



2l/6r(|) ^ "^ ' 32/3r(|)(l/ + 02a)l/3' 

hence 

^ V2^ ^ V2v^/^ / _ 2^ 

^^ 3(z/ + 6i2a)i/3 3 V 3i/. 

Thus, using this and setting C3 = C2 — Ci, we can write 

(36) y{z) = 21/3^.1/6 v4z(-2i/3z)- 

^^ Ai{-2''^z) + C3 Vi Ji/3(V2 C) + e2i?(z) := 2'/^v''^ Ai{-2''\) + 7^. 
Now Ai(— 21/^7) — and by Lemma [T2l we have 



w ' 3 20 3vi^ 

yielding 

10 ,^ /- ^7T, 40 

C3 < -^ (2 + VI + iv-^/^) < 



3v^ 3a/j7 
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On the other hand, 

giving C3 > -7/v^. 
Inequahty (^01 ) gives 



7-27/12(4^ + 301/3)1/4 ' 

and using ^ and p^ we have 

VI|Ji/3(V2C)| < -M^^^ , 2V4vlzV4} • 

Combining these estimates with psp after some straightforward calculations one 
finds 

23max{M^ 
' ' "^ 2v^ 

and 1331 follows. D 

6. Sharper asymptotics 

The classical asymptotic given by (l?fl) does not makes much sense for x — 0{fj.) 
when the main term and the error are of the same order. Here using formula ([7]) 
we derive a different asymptotic expression with much smaller error term. It also 
leads to very sharp approximation of the Airy function Ai{—x). 

We'll need a few lemmas given in 18^. 

Lemma 14. Let f{x) satisfy the differential equation 

r'ix) + b\x)f{x)^o, 

where b{x) > and b"{x) exists on an interval T. 

Let g{x) = y/b(x)f{x), 

then for x £ I, provided the integral exists, 

3b'^{t)-2b{t)b"{t) 
463 (t) 

where B{x) = J b(t)dt, a £ I is arbitrary and \9\ < 1. 

Proof. Observe that g{x) satisfies the equation 

(38) g"-jg' + g{b'-e)^0, e ^ e{x) ^ — . 

The solution of the corresponding homogeneous equation 

is 

go — ci sinS(x) + C2 cosB{x). 

Solving formally p8p as a nonhomogeneous equation with the right hand side 
e{x)g(x) we get 

g{x) = go(x) + j ''-^^ ^ e{t)g{t)dt = 



(37) .g(a;) = ci cosS(x) + C2 sinS(a;) 



9{x) 



dt. 



20 



goix) + 9 
go{x) + 
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sm{B{x)-B{t)) 



m 



it)9{t) 



3&'2(i) - 2b{t)b"{t) 



4fe3(i) 



'W 



dt = 
dt. 



The normal form of differential equation ^ is 

/" + (1 - ^^-^)/ = 0, f^yGMx). 



Thus, for X > 1 /max{0, v"^ — j} we have 



6(a;) = 
\/x^+7i + y/]^ In 



2-2 _ j^2 _|_ 1. 



B(x) = 



^+^/|J.+x^ 



and 



y/x^~~jL + y/Jl arcsin ^ , 



1^1 < 



X > 



Theorem 15. For |zy| < i anrf a; > 0, 



(39) J^{x) ^J-ix'+ ti)~''^ cos {B{x) - uj^) + 

i^or \v\ > i and x > fi, 



12-KX (x2 +/i)3/2 



13/i 



(40) J,(x) = J-(a;2-/,)-V4eos(S(a;)-c.,)+L 

V TT 12V27r (a;2 - iifl'^ 

Proof. Since | Ji.(a;)| < J^ for |z^| < i , by dS?]) we have 

6z^ + ^ 



g(a;) = .90 + ^ 



M 



78^7, z3/2 (z2 + /i)9/4 

Comparing this with the standard asymptotic 

TT 

for large x one finds 

'2 



dz = go 



(41) 



/(x) = ^/xJ^{x) 



cos(a;-a;i.) +0(x"^/^), 



Cl 



smw^, C2 



COSWj^, 



yielding (jM)) . 

Similarly, for j/ > i and x > /i, using p7|) instead of (fT3|) . we obtain 



D 



9{x) =.90- 



8^ 7, (z2 - ^)5/2 



6z — /i 1 

^7:r dz = .90 + ^ 



Stt 



3x2 ^_ 2^ arcsin -^ 



3(x2 - ^)3/2 ^ 
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It is easy to check that the function 



3a;2 + 2^ arcsin ^ 



3(2;2 - ^)3/2 ^ 

is decreasing and positive. Therefore, by arcsin— > — , x > 0, we get 
3x2 ^ 2fi arcsin ^ 3x^ + 2^ 1 _ 



3(x2 - /i)3/2 ^ - 3(a;2 _ ^)3/2 a; 



13Ai (x + 2^x2 - A* )( Va;^ - fi-xf 13Ai 



6(x2 - /i)3/2 2x(x2 - ^)3/2 6(x2 - /i)3/2 ' 

and 

0(X) = On + t^ z= — . 

^^ ' ^ 12^/2^ (x2-/i)3/2 

Comparing this with the asymptotic for large x one finds 

ci — \ — smwy, Ci = -i/— coswy, 

V TT V TT 

and dlOl) foUows. D 



Remark 2. T/ie approximation given by ^JW *s strong enough to be matched with 
the solution in the transition region, as for x — v + v^'^ z the error term is of 
order 0{v~^'^ z~'''^). This allows one to find the constants of integrations ci, C2 in 
Theorem \13\ avoiding delicate claims like Lemma \JM We omit the details. 

Similarly to Corollary [T] we obtain 

Corollary 2. 

(42) Ai{-x) = 



v^ (16x3 + 5)1/4 0Fxi/4(16x3 + 5)3/2 ' 

provided x > 0. 

The argument of the cosine in (|5^ and (|42l) (but not in (HH])) can be simplified 
at the cost of a weaker numerical constant at the error term. Namely, it is easy to 
verify the following elementary inequalities: 

2 

y/x^ + ^i + ^ In ^ "", J =^"^ + ^'^' ^>0, 

y/JJ+y/fi + X^ 2x 24x3 



Vl6x3 + 5 VS , Vl6x3 + 5 + V5 _2 3/2 _ _5_ -3/2 n2^^ -9/2 
6 6 4x3/2 3 48 9216 

Since | cos(x + e) — cosx| < e, we obtain 

^ ' ' Vtt (x2+^)i/4 24V2^x3(x2 + /i)i/4' '11-2 

2^/icos(-x3/2 - A2:-3/2 _ IL\ K 
(44) Ai{~x) = ^ '^3 4^^ 1^ ' « ^ 



0F(16x3 + 5)1/4 9V7Fx4(16x3 + 5)1/4 



22 I. KRASIKOV 

In particular, the last formula yields rather sharp approximations for the zeros 
of Airy (e.g. already for the first zero the error is less than 0.00122), and, in view 
of the inequality dSH]) . the Bessel function. 

Theorem 16. 

where as is s*'' positive zero of Ai{—x) and m = (12s — 3)7r. 

Proof. Elementary arguments show that |sina;| < e implies sin(x + ^^) — for 
some 6, \9\ < 1. Therefore, Ai{—x) = means 

I ^'^' - Js ^"'^' - T + ^ S ^"'^' = "'' ' = °' ^' ^' - 

Since ai — 2.33..., we may assume x > 2. Then 

36 



Aa;-3/2 > ^3,-9/2 



hence 

Z 0/9 OTT 

_a;J/- > _+7rs, 
and we obtain the following estimate 

^^-9/2^ 640 



36 2437r2(4s + 3)3 ■ 

Thus, for (s + 1)*'* zero this gives the equation 

2^/9 5 o/n Stt ^ 6407r 37r 

- x''/'^ x"-^'-^ = h TTs + 61 :^ h ITS + e, 

3 48 4 9m3 4 ' 

with the relevant solution 

X = 16"^/^ fm + 12e + ^(m + 12e)2 +40] 
After some calculations one gets 

a; = 16"^/^ (m + Vm^ + 4o) (l + e^=^^] = 
V / V Vm2 + 40 / 

>2/3 „ 12807r 



16"^/^ [771+ ^^2+40] +( 



9TO3(m2+ 40)1/6 ■ 
This completes the proof. D 



Formula p5|) can be simplified at the cost of slightly weaker numerical constant. 
Namely, as one can checks 



< -(^2 + 20)^/'"^ - 16^2/^ (m + Vm^ + 40 ] 



2/3 25 

< 



3to3(to2+ 40)1/6 
yielding 



4^ m3(m2 + 40)1/6 ' 

Finally, comparing the numerical values of the zeros of Ai(—x) with (P5|) leads to 
the following conjecture: 

as < 16-2/3 (m + v/m2 + 4o] . 
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SOME ASYMPTOTICS FOR THE BESSEL FUNCTIONS WITH 
AN EXPLICIT ERROR TERM 
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Abstract. We show how one can obtain an asymptotic expression for some 
special functions satisfying a second order differential equation with a very 
explicit error term starting from appropriate upper bounds. We will work out 
the details for the Bessel function Jy{x) and the Airy function Ai{x) and find 
a sharp approximation for their zeros. We also answer the question raised by 
Olenko by showing that 



u^- 


1 


< sup x^/^ 


Mx) - 


J — cos (x - 


TTl/ 




4 


x>0 




\ TTX V 


•2. 



< C2 



i^ > — ^ , for some explicit numerical constants c\ and C2. 

Keywords; Bessel function, Airy function, asymptotic, error term, zeros 
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1. Introduction 

All basic formulas and asymptotic expressions for special functions we use with- 
out references can be found in [12]. To write down error terms in a compact form 
we will use 0, 0i, 02, ■•■, to denote quantities with the absolute value not exceeding 
one. 

In most of the cases error terms of asymptotics of special functions are either not 
known or, at best, valid for a rather restricted rang of parameters. The following 
is a typical example of that kind (see e.g. [12] Ch. 10]). 

The Bessel function J,\x) is defined by the series 



(1) 



U-) = (f ) E(-i)^- 



(a;V4p 



i=o 



j!r(j 



1) 



and is a solution of the following ODE: 

(2) x^J'lix) + xJl{x) + (x? - v^)J^[x) 



Theorem 1. Suppose that i' > 0, x > 0, lUi, 



2 



then 
(3) 



> maxf ,1), ii> max( — 



— Jy{x) = cos (x - Ui,) 



0. 

and let 
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sin(x - ^.) I 2^ ..+1 + ^^ ^^^^ 



, i=0 



a;^ 



2 2;2^2 + l 



where 



ai{v) 



2HI 



The assumption i^ > is not really restrictive and can be surmount by, say, 
applying the three term recurrence for J^. However if p is large or depends on x 
estimating the error term in ^ seems at least as difficult as the original task. 

In this paper we show that there is a simple way to circumvent this problem 
and to find an explicit expression for error terms which is also uniform in the 
parameters, provided one has an a priory upper bound on the absolute value of 
the considered function. In turn, in many cases such a bound may be obtained 
by using so-called Sonin's function. For Bessel and Airy functions, as well as for 
Herniite polynomials (see [2]), this can be done in a quite rutin way. For Jacobi 
and Laguerre polynomials it is a much more involved problem and the result is 
known only for oscillatory and transition regions [6], [7], [8]. It is worth noticing that 
despite the fact that it is rather a technical problem and we do have appropriate 
tools to tackle it (see e.g. Lemmas [H and S] below), one still needs a good deal 
of calculations to extend the bounds to monotonicity region. Thus, although the 
underlying idea of the method we use here is quit simple and can be applied to 
many special functions satisfying a second order differential equation, it is not 
utterly straightforward to work out the details. In this paper we will consider 
the Bessel function Jv{x) as an important example to illustrate this approach. 
We provide asymptotic expressions with an explicit error term for oscillatory and 
transition regions and also give some new estimates in the monotonicity region. As 
a corollary we derive a surprisingly accurate approximation for the Airy function 
Ai[—x)^ X > 0, and obtain an approximation with an explicit error term for its 
positive zeros, which, due to known inequalities, yields an approximation of the 
zeros of the Bessel function. 

It is worth noticing that the obtained estimates are, in a sense, best possible. 
In particular we will answer a question raised by Olenko [llj by showing that for 



v> --. 



(4) c. 



< sup a; 

x>0 



3/2 



Jty{x) 



2 / nv n 

cos [x 

TTX \ 2 4 



< C2 



for some explicit numerical constants ci and 02- 

The paper is organized as follows. In the next section we describe an idea of the 
method. In section [3] we establish some upper bounds on Bessel and Airy functions 
we need in the sequel. Our main tool here will be Sonin's function. In section |4] we 
consider the error term of the standard asymptotic 



M.)-^.o,{.-'i-'-). 



and prove ((4]), thus answering Olenko's question. The error term for asymptotic in 
the transition region will be derived in section [5j Finally, in section [51 using the 
approach of section[21 we establish a different sharper approximation for Bessel and 
Airy functions and their zeros. 
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2. Preliminaries 



The following approach was shortly described in [8] . We want to find an approx- 
imation of a solution of the differential equation 

(5) f" + b^{x)f{x)^0, 

in terms of some standard function F{x)^ which also satisfies a second order ODE 

Vi{F)^P2{x)F"+pi{x)F'+po{x)F^0. 
In fact, in what follows we choose F to be just cos (pix) with an appropriate function 

A possible approach to this problem is to seek for a multiplier function z{x) such 
that the differential operator 

2^2(5) = q2{x)g" + qi{x)g' + qa{x)g 

for g — g{x) — z{x)f{x) is in some sense close to I?i. For example, in the WKB-type 
approximation one chooses g{x) = \/b{x) f{x), yielding 

h' W^ ~ 9hh" 

(6) 9"~'-g' + b'gil + e{x))^0, e(x) ^ J'' . 

If e is small we can expect that g{x) is close to the solution of the equation 

h' 

" / , 1,2 n 

9o --^90 + 90 = 0, 

which is just go = McosB{x), where B{x) = J b{x)dx. 

Assume now that we have an a priori bound 15(2:) | < C. Then we can readily 
estimate the error term Ig — gol by solving ([6]) as an inhomogeneous equation. 



(7) g{x)^go-J e{t)b{t)sm{Bix)-B{t))g{t)dt, 

thus obtaining 

g = go+eC I e{t)b{t)dt. 



To derive an upper bound on |(7(a;)| we consider Sonin's function 

^'2 452 

S{x) = S{g; x)^g^ + ^.^^—^ = g- + 4^4 _ 256" + 35'^ ^"' 



then 



, _ 8b{6b'^ - 6bb'b" + b^b'") ,2 

^ y^) - /AU4 , OU2 nuuf\2 9 ■ 



(464 _,_ 3^/2 _ 255") 2 

Thus, if 46* - 266" + 36'^ > and 66'^ - 666'6" + 6^6"' > then S" > 0, and we 
obtain .g^(x) < 5 (00). 

Moreover, one can also get an upper bound on S in the following way: 

__ 6(464 + 36'^ -266") 

2(66'3 - 666'6" + 6'26"') ^ " ' 

that is 

2(66^3 _ 6j,yy/ ^ y2y//) ^ ^ ^ ^4 

■^ /'^ - (./^M , 01./?. o^,^,//^ ~ ZJZ '""" 



6(464 + 36'2 - 266") dx 464 + 36'2 - 266" 
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provided the last expression is nonnegative. Integrating, we find 



Siy).^ + m 


b"ix) 
263(2) 


1 + e{x) 


^(-)^l+l^ 


b"(y) 

263(a) 


1 + e(y) 



This shows that the envelop of g'^{x) given by S{x) is almost a constant as far as 
e(x) = o(l). 

In practically important examples the situation is somewhat more subtle as 
the coefficient b{x) may vanish. For instance for the Bessel function b{x) — 
x~^ ^x^ — i/"^ + 1/4 , and Sonin's function does not provide any information for 
the monotonicity region < a; < -^/i/^ — 1/4. Thus, one needs some supplemen- 
tary estimates to extend the bounds on \g{x)\ to this interval. Let us notice that 
although the behaviour of the solutions of (O looks less complicate in the mono- 
tonicity region, it, probably, allows only a piecewise approximation in reasonably 
simple elementary functions. 

Another rather technical problem is how to find the constants of integration in 
go- Here one either has to know, at least approximately, the value of g(x) at some 
points, e.g. at infinity, or to be able to match asymptotics in the oscillatory and 
transition regions. 

It is worth noticing that an asymptotic with an explicit error term in the transi- 
tion region, that is around a zero of b{x), can be obtained quite directly, provided 
we know a bound on |g(a;)|. Indeed, let b{a) = 0, and let d = ^ b'^{x)\ _ , then 
in a vicinity of a we can write 

b''{a + d-^/h) = <f^h + S{t), 

where 5{t) is small. The function y{t) = f{a + d^^/'^t) satisfies an Airy type ODE 

Solving it as an inhomogeneous equation one gets an explicit error term as above, 
yet facing again the problem of fixing integration constants. 

3. Upper bounds 

For the Bessel function it will be convenient to introduce the parameter /i = 
I i^^ — I; I . We also use the above notation w^ = ^ -I- ^ . 

In this section we establish some upper bounds we need in the sequel. A simplest 
inequality of this type [TH] states that for x real 

(8) JAx)< 



2Tfi/-|-l) 



For our purposes we need much more accurate estimates. To bound the Bessel 
function in the monotonicity region we will apply the following inequality given in 
[5] . We sketch a proof for self-completeness. 

Lemma 2. Let Ju{x) = x"" J,y(a;), v > — ^ , then for < a; < i^ -I- ^ , 



(9) JM. > V(2^+l)^"4x^-2z.-l > 2x 



Jt^ix) ~ 2x - 2iy + l 
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Proof. The function Jv(x) is an entire function with only real zeros satisfying the 
Laguerre inequality J'J' — JuJU > 0. Substituting here J" from the differential 
equation 

xX' + {2v + l)Jl + xj^ = 0, 
and dividing by Jy jx we obtain 

xt^[x) + {2v + \)t{x) + a; > 0, 
where t{x) = Jlj Jv Hence, 



, , , , / ^{2v+\f-'^x'^ + 2v + \ 2x 



2.T y/{2iy+iy -4.t2 + 2z/ + 1 



Since t{x) — ?► when a; — >■ 0+, whereas — ^^^- g^ > — oo, we conclude 

that 

2x 2x 

t(x) > = > 



v/(2z/+l)2-4a;2 + 2j/+l " 2zy + 1 ' 
We need the following very accurate estimate giving the value of Jy{i') 

where a = 0.09434980..., [J. 

The following theorem improves the inequality 

J^itv) <Jy{iy)fe'^^-'^'', iy>0, 0<t<l. 

given in |13j . For large v it is also stronger than the classical inequality 

/ (r) < — p-^V4(^+i) 



D 



PI p. 16], provided t > y/\n 16-2 w 0.88. 
Theorem 3. For x = tv, < i < 1 and v > 0, 

(11) JuM < Miy)^ exp ,\ / < ,,,,^,2, „^,,^ exp 



2z^ + i y 32/3r(|)t/''+i/3 ^\^2i/ + i 

Proof. By the previous lemma we have 



Jy{x) - J^ 2v + l 2v^\ 

hence 

Mty) < J.(j^)exp ; , / 
V 2i/+ 1 

which together with ([T0| yields the required result. D 

Remark 1. The function J^,j{x) = x^'^J^^x) of Lemma\^ belongs to so-called 
Pdly a- Laguerre class and satisfies the infinite series of inequalities: 

2m (2m\ 

(12) L„,(J,) = Y.{-ir+l^jU)j(^rn-,) > 0^ 
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where Li is the usual Laguerre inequality J'J' — JvJ'J > 0, (see e.g. |14j . [15j ). 
Using Lm,{J^ii) > for m > 1 leads to much more precise yet more complicated 
bounds on Jlj Jy and consequently on J^. Alternatively, one can use the inequality 
Li{J^u + ^J'l) > 0, A G M, then optimizing in A. It is worth noticing that both 
methods give an inequality .similar to (0) but in the opposite direction. Thus, one 
can use the known value o/j7iy(0) instead of Jjy{i>). 

Our main tool for bounding solutions of the second order differential equations 
will be Sonin's function. In particular, it was used by Szego to prove that for 

(13) \Mx)\ < 

V TTX 

Although he did not state this explicitly, his proof of Theorem 7.31.2. immediately 
implies 

(14) \YAx)\<\f^- 

V TTX 

His arguments go as follows: let y be a solution of the Bessel differential equation 

x'^y" + xy' + (x^ - iy'^)y = 0, 
then for |i^| < | Sonin's function 

(15) '^^^-^y' + J^.ii-^^y^ 

is increasing and inequalities (jl3|) and (jl4[) follow by calculating S{oo) from known 
asymptotics of J^ and Y,y, whereas for i/ > ^ it is decreasing for x > ^ffl, and does 
not lead, at least directly, to any explicit inequality. 

It turns out that for v > 1/2 it is more natural to deal with the function 

(16) H.(x) = |a;2-/i|i/V,(a;), 

rather than y/xJy{x). Here we will refine an inequality for the Bessel function 
obtained in [2]. First we need a bound on location of the leftmost maximum of 
H,{x). 

Lemma 4. The first positive maximum ofT-Ly{x), v > \ , is attained at a point ^ 
satisfying 



^>v^l- (2:/)-2/3 



Proof. Since obviously < ^ < /i, we can restrict ourselves to the interval (0,/i) 
and write down 

n,{x)^x''{pi-x^f'^Mx), 

where as above J,y{x) = x^" Jy{x). Then 

= n'M) = 2{f-eY/^ ^^^^ " ^')(^^^^^ + ''^ " ^'^ ^"^^^^ 

where t{x) — Jl{x) / Ji,{x) . Hence 



m 



(4j/2 - 1 _ 4^2)^ 
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and comparing this with ()11|) we obtain the inequahty 



4i/2 - 1 - 4^2 - 2 ' 

Simphfying we get 

16f-Ai2iy+l){6iy-l)^'^+{2iy~l){2v+lfi6v+l)S.^+iy{:y+l){Av'^-lf ;== p(0 > 0. 

Observe that for v > 5/3 this polynomial has the only positive zero ^07 Indeed, the 
discriminant of p(f) in ^, up to an irrelevant numerical factor, is 

v{v + l)(2zy - l)6(2zy + l)i°(108i/2 - 172zy - 5)^. 

Thus the number of positive zeros does not change for 108i^^ — 172i^ — 5 > 0, 
in particular for v > 5/3 . For v = 5/2 we obtain the following test equation 
^6 _ 21^4 ^ ;^44^2 _ 325 = 0, with the only positive zero S, w 2.14. Finally, 

p{,y) = u{'iv'^ - 2v ~ I) > 0, 

and using the substitution n = r^/2, wc find 



p{v^Jl - (2z/)-2/3 ) = p(v/r6-r4/2) = -r^r^ - lf{2r^ + 4r^ -r + 2)<0, 
hence 



e>I^Y^l-(2l/)-2/3. 

D 



Theorem 5. Let v > ^ , then for x > 0, 

(17) \x' - f,\'/VA^)\ < 

and the constant J ^ is best possible. 

Proof. For x ~ ^/Jl the result is trivial. Otherwise we shall consider three cases. 
Case 1: X > JJL. The function 

H.(a;) = (a;2-/i)i/4j,(x), 
as easy to check, satisfies the differential equation 

Ki^ - ^^—, HUx) + 4(-^ -/;);+ (6^;- A-)^ H.(x) . 0. 



Consider Sonin's function 



4x^(a;2 - /i)2 



^^") ^ ^' (") + 4(.2 :;)r+ (eS _ ^)^ ^'' (^)' 



then Ul{x) < S{x) ioi x > ^ > 0. 
One finds 



2Afix^{x^ - lj.){Ax^ +^) 



a, 



s'{x) = — ^— ^ !^ ^n'tix) > 0, 

(4(a;2 - ^)3 + (6x2 - /i)^)2 ..w - , 
hence 

\HUx)\ </\im S{x). 

Using 

w. . Ji.-i(a;) - J,y+i(a:) 
J AX) = ^ 
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and the asymptotic formula 

Jyi.^) ~ \ — cos(x ), 

V Tra; 4 

after some calculations one finds 

Since 'Hl{x) — S{x) at all local maxima the constant \/ ^ is sharp. 
Case 2: < X < y/Jl, 5 < J^ < 3. By dH]) and ly > ^ we have 



nAx) = i,-xyf^Mx)<^jy^^ 



The maximum of the last expression is attained for x = ^/v^ — v/2 , yielding 

For i^ > 5 the function T{v + 1) is increasing and it is easy to check that f{v) is 
also increasing. Hence for 2 — '^ — ^ — ^^ 0,7^ 6, we have the following estimate 

H.(x)<-^. 
r(a+ 1) 

This yields 

35/4 35/4 r^ I 

4V2r(|) 2^/2^ V TT ' 2 

2 72 [2 

H.(.)<j^<^-, 2.6<.<3. 

Case ■?.■ < x < ^, i^ > 3. Inequality (|TT|) yields 



^-^"^< 3V3r(|K+i/3 ^"PV2.+ 1 



Let r = (i^/2)^/3^ ^y Lemma [4] we can set 



X = vJl - {2v)^'^/^z = — \Jr'^ - z , r "^ < z < 1. 
This gives 

n.{x) <a{x- ^f [zr^ - 1)VV- exp (^) := A/(z), 
where A — 3^3 ,^, . We find 



/'(z) _ (l + r3+r'^-2zr'^-zV)r3 
/(z) ~ 4(l + r2)(r2-z)(zr4-l) ' 
where 

l + r'^ + r^ -2zr^ - z^r^ > (r^ - l)(l + r)(r3 - 2r'^ + r - 1) > 0, 
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for V > 19/7. Hence f{x) is increasing and 

nAx) < Am ^a(i-^J^' (r^ - l)i/V-i exp (^-^) < 

Ae-^/^(l-.-)V^exp(^)<A</f. 
This completes the proof. D 

A similar but more complicated result can be given for J^{x). 
Theorem 6. Let v > \ and x>v+ J1J3 v^l'^ , then 

x{A{x^-i.^f-ix^-lQx^v^ + ^^Y'\ ^,^ ^, ^ 2 

X V V ^ 

Proof. Let 

x%l)'^/'^(x) 

z{x) = — ^ J^ X , 

X'^ — v^ 

where 

V'(x) = 4(x^ - v^f ~ Sx** - 10x^2 + i^"*. 

First notice that '^(x) > for i/ > ^ and a; > z^ + -^i/T" i^^^'^. Indeed, the substitu- 
tions 

(19) :, = y + . + ^^.V3^ . = (2-1/3 + n)^ 

transforms V' into a polynomial in n and y with nonegative coefficients. Consider 
Sonyn's function 

Si^x) = -Ax)^^^^^^^^.'^ix). 

and its derivative 

nx) ^ ^-^(-^ - ^:)>(^)^(^) ,,.(,), 



where 



g(x) = 16(a;2 - 1.2)9 _ 4(3.2 _ ^,2)6(2^4 _ 24iy2a;2 _ 9^.4)^ 

{x^ - v^f{v^ + ^l^v^x^ + \mv^x^ + 45x*)+ 
3a;2(2iyi° + 2^v^x^ - lUiy^x^ - 66j/'*a;^ - 42iy^x^ - Sx^^), 



and 



P{x) = 

\%v^(x?- - v^f{^x^ + v^) + 8(a:2 - i.2)-''(5i.^ + 119i^V + SISj^V + \\hv^x^ + 6i/®)+ 

17!^i2 + ^\%v^^x^ + 695(^*2;'' + 492iy^x^ + 2831i^''a;^ - h^v^x^^ + ^x^'^- 

2v^{v^ + Sli^^x^ + 'iiv^x'^ + ITTi/^a;® - 18a;^). 

Applying ([T51) to P and Q we obtain polynomials with nonnegative coefficients. 

Hence, z'^{x) < S{x) for z^ > i and x > v+ 2^/3"'" '^^Z'^. Finally, using the asymptotics 

Jy{x) ~ W — cos(x -w,,), J'(a;) ^ ~\ — sin(a; -Wi,), 

V TTX V TTX 
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we obtain 

z^(x) < lini S{x) ~ — , 

a:— 7-00 TT 

and the result follows. D 

We'll need one more statement of this type for the Airy function 

M-X) = ^ (J-1/3(C) + Jl/3(C)), 



where ( — ^ — . 
Lemma 7. 

(20) {x + c)^^^Ai{-x)<—, x>-c, 

where c = 15^'^ ■ 2"^'"^. Moreover , the values of all local maxima of the function 
{x + cYi'^Ai{—x) are restricted to the interval {—i= , jg). 

Proof. Consider the function 

f{x) = {x + cy^^Ai{-x), x>-c, 
which satisfies the following differential equation: 

^''(^)-2(^^'(^)+(^+16(^)^(^)-«- 
The corresponding Sonin's function is 

S{x) = /2 + , ^ , 

^+ 16(c+a;)^ 

and 

, , , 256 cx(x + c)(x + 2c) 

b (x) ~ ^— . 

{16x{x + cy+bf 

Hence, a; = is the only maximum of S{x) for x > — c and 

fix)<f{0<^, 

where ^ = 1.12879717..., corresponds to the first maximum of f{x). Using the 
asymptotic M(— x) ^ n^^/^x^^/^, where M{x) is defined by Ai{x) — M{x) sin(/)(x), 
one finds 

lim f{x)=TT-^^^, 

a:— >oo 

and the result follows. D 

In the transition region the following inequality, which may be of independent 
interest, will be useful: 

Theorem 8. For xi,X2 > and < ly < ^, 



(21) y/x^ \J^i,{xi)J^{x2) - J^,y{x2)Jiy{xi)\ < — siuTTJ/. 

TT 
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Proof. The function 

F = F{xi,X2) = y/x^ {J-u{,Xi)Ji,{x2) - J-v{x2)Jv{xi)) 

satisfies the Bessel differential equations 

(22) A_F + 5(a;,)F = 0, i = f,2; 



where 






For a; > we consider the following majorant of F given by Sonin's function 



(^F\^ 

F^XI,X2)<S{X1,X2)^F' 



b{x2) 
By (122]) we have 



_c _ dX2 ^ ^^ ——F 

dx2 b'^{x2) \dx2 

Since 

8^2'^^'^ = ^^^' 
for 1/ < ^ , the Sonin's function is increases in 2:2. Using the asymtotics 

J,(:,) = ./X cos(x - (^!:±1)!L ) + o{x-^/'), 

\l TTX 4 

•^^ a^ = - V — sma: - ^ -^ + o x-1/2 , 

V TTX 4 

and inequalities (|13l) . ((Ti)) . on taking the limit we obtain 

F'^(x\,X2) <S{x\,X2) < lini S{x\,X2) = 

3:2—^00 

(j2(xi)-2cos7riyJ^(a;i)J_^(xi) + J2^(xi)) = 

2xi sin "Kv , ^^ , , ,^9, ,^ 4 sin -kv 
^— (J2(xi) + Y^{x^)) < ^, . 

This completes the proof. D 

To get better numerical constants the following inequalities will be useful. 
Lemma 9. For a: > 0, 

r Isintl , 2 
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Proof. We have 



sm^t , 1 1 /■°° cos2i . 1 



dt= / TTT dt < —- 



fi-)--l TT^dt. 



/o {t + x)2 2x 2 Jo (t + a;)2 2a; 

because 

f°° cos2t , ^ f^^^^ cos2t , „ 
/ 7 ^ dt = > / 7 r^ dt > 

io (t+xr f^j^ {t+x)^ 

is an alternating sum with decreasing terms. This proves 
Now we prove p4)) . Let 

I sini| 

(t + xy 
We have 

f^ xt , f°° xdt , , x + l 2 ,3, 

Next, we show that 

2 
hm f{x) = — . 

a:— foo 7r 

Using I sin(a: + 7rfc)| = sinx, < x < § , we obtain 

T(k+1) I • ,1 oo „7r • , 

dt = 






(i + a; + ttA:) 



2 



^X ''"'^„(^T7T^^'^'^^y„ V'(^)smtdt, 

where ipi^) — T'{x)/T{x) is the digamma function. Since ip'{x) = x~^ + 0{x^^) 
for a; — > oo this yields 

X f^ sint ,. „, _i X [^ smt _-^ 2 _i 



f{x) = - / dt + 0{x-') = - dt + 0{x-') = - + Oix-'). 

n Jq t + X -K Jo X TT 

Now it is enough to show that f{x) is increasing for x > j . Using the inequalities 

1 1 ,// ^ 1 1 1 



1 + X X^ 1 + X x^ {x + 1)2 



we have 



^7' (a;) = ^ (x r ^'{l^^) sintdt, 
dx \ Jq n 

V''( -)smtdt+- %}:"{ ^)sintdi = 

"■ '^ Jo '^ 

^ t + X X 1"^ t^-x 

V''( ) sinidi / ^' { )costdt> 

^ 7^ Jo I" 

2 

sin t dt 



V^ + ^ + '^ (t + x) 

' ' X TTX TTX , 

COS t dt 



t + X + TT {t + xY {t + X + Tl) 

X TTX 
h 7" TTT 1 COS tdt > 

jr/2 \'t + X + TT (t + xy 
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(j + l)7r/4 

sin t dt 

JTr/i 



costdt 



where 



2Q{x) 

7r-"P(2;) = 



f' 



1146882/" + 4096(259 - 2^2 )y^" + 1024(4216 - 77\/2 )y^ + 1536(6641 - 244\/2 ) 
64(242859- 16954^2)/ +400(40011 - 5002^2)/ +4(2820158- 591757V2)y5+ 
4(1332737- 443447 V2)2/'' + (1564982 - 808187V2)y3 + 5(44796 - 40571^2 )2/2- 

3(1140 + 7109V2)y- 3780, 
and 
Q(x) = y(y + l)2(2y+l)2(4y + l)2(4y + 3)2(y + 2)(2y + 3)(4y + 5)(4y + 7), y = x/tt. 

It is easy to check that the polynomial P{x + I ) has only positive coefficients. 
Hence f'{x) > for a; > 2/3, and f{x) < liiRx^oa f {x) — - . This competes the 
proof. D 

4. Oscillatory Region 
Having at hand an upper bound on | Jiy(a;)| one can estimate the difference 

/ TTX 

(25) r{x) = J — J^{x) - cos{x - uj^), 

in a rather elementary way. Notice that r{x) satisfies the following differential 
equation 

with the general solution of the form 

[^ f'^ ?\ f°° sm(t - x) ^ , , , 

(26) r(a;) = CiCosx + C2Sinx + W- [j^'^ I / — ttt^ J,y[t)dt. 

Now one has only to estimate the integral and to notice that as far as it is o(l), we 
have ci = C2 = by an obvious limiting argument. 

In |llj Olenko proved the inequalities which for v > can be written as 



ciz/^/^ < supx^/2 

x>0 



Ju(x) - \l — cos(a; - uj^) 
" nx 



< C2V 



13/6 



with explicit constants ci,C2, and raised the question what is the best possible 
exponent a oi v h\ these inequalities. It turns out that that the answer a = 2 can 
be readily extracted from (|26p , since starting with a reasonably sharp approximation 
to the Bessel function one can iterate it getting more and more accurate yet more 
complicate approximations. 
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Theorem 10. 

(27) 
where 



Jv{x) = \ — cos(x — Wy) + Oc^x 



-3/2 



f {^" > 



a;>0, 



k' < 



2 



X> y/JI, V> 



|, 0<x<^, z/>i. 

Moreover, up to the numerical factor c the error term in {21^ is sharp. In particula 
c cannot be taken less than l/yjzn . 

Proof. To estimate the integral in ([25)1 for |;/| < ^ we apply (|24p yielding 



I^{x) 



sin(t — a;) 



J^(t)di 



< 



|sin(t-a;)| , 2 



i2 



Tree 



Thus 



2/i 



r(a;) = ci cos x + C2 sin x + i 

t:x 

For 1/ > i and a; > y/JI we use ([TT]) and the inequality arcsina; < ^ . This gives 
I sin(i — x)\ 



l^x) < 



t3/2(t2_^)l/4 



dt < 



sin (i- x)dt 
t^ 



dt 



X tv^ 



' arcsm 



Vm 



< 



y IxJJL 2x 

Hence in this case 

r(x) — ci cos X + C2 sin x + 9 
Similarly, for j/ > ^ and < a; < ,/Jl , 



2x 



I^{x) < 



x/M 



I sin(t — a;) I 

t3/2(i2_^)l/4 



dt 



\/m 



I sin(t — a:) I 

, t3/2(/i-t2)l/' 



■dt < 



1 /¥ 
2Vm 



v^ |sin(t-a;)| 1 /¥ 

., i3/2(^_i2)l/4«^- 2 Y/i' 



t2 



x/m 



dt 



Wli - 1^ 



1 /F 
2Vm 



2x 



2^ 



/iX X 'ZX Xy/JI 

An elementary investigation shows that the maximum of the function 

x(VM-.)ln^ 

X 

is attained for x = d^/JI , d — 0.314711..., and does not exceed 2///5. Hence, 

I.(x) < ^^ < 
and 



r{x) — C\ cos X + C2 sin x + 



7x ' 



7x 
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Finally, since for a; — >■ oo, 



J^{x) = \ — cos{x ~ Lu^) + 0{x '^/^), 

V TTX 

that is lima;^oo r(x) — 0, we conclude that ci — C2 — and (|27|) follows. 

Let us show now that up to the numerical factor c the error term in (j27p is sharp. 
By (HH) and ^ we have 



n^{x) =x^/^ 



Jv{x) - \ — cos(a; - u^) 
nx 



lix 



sm(t — x) 



t3/2 



J^{t)dt 



— iix 



sin(i — x) cos(i — uj^) 



t^ 



dt 



-6cfix 



sin(t — x) 



P 



dt 



:= \l - fj.x\Ii\+l2- 



Here 



I2 — Ocfix 



dt Ocfj, 



To bound /i we introduce two auxiliary functions / and g (see [HI Ch. 6]), defined 

by 

TT 

Si{x) = — — f{x) COS a; — g{x) sin a:, 



Ci{x) — f{x) sin a; — ^(a;) cos a;, 



with the asymptotics 



/(a;) = - + 0(a; ^), g(a;) = ^ + 0(a; ''), x -^ 00. 



Calculations yield 

/i = 5*1(22) sin(x + cji,) + Ci{2z) cos(a; + uj„) - 

f{2z) sin(2i + a: — cui,) — g{2z) cos(2t + a; — uj^) 
sin{LL>i, — x) 



sintcos(2: — Wj/) 



sinicos(z — lUi,) 



t=o 



2x 



+ 0{x ^), z^x + t. 



Hence 



M 



K^{x) — — \sin{LJiy — x)\ + 0{x ^), 



/27r 



and the result follows. D 

Applying (j27p to the Airy function yields the following approximation: 
Corollary 1. For x > 0, 



(28) 



Ai{-x) 



cos(C - f ) 



/^a;i/4 6\/3 7r3/2a;7/4 ' 
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5. Transition Region 

Our estimates in the transition region are based on the following simple obser- 
vation. 

Lemma 11. Let for a real constant q, 

y"{x) + q xy{x) — u{x). 

Then for a; > 0, 

(29) y(x) = V^ j^ci J_i/3(^i^— ) + C2 Ji/3(^^) j + 

9q-^x-^/^ r \u{t)\t-^'^dt, 
Jo 
provided the integral exists. 

Proof. Let yi and 2/2 be two linearly independent solutions of the homogeneous 
equation 

y"{x)+q^xy{x)=Q, 

and let 

rj, ,N ^ Vi{i)V2{x) -~yi{x)y2{t) 

^'' yiiMit) - y'iit)y2{t) - 

Then 

y{x) = ciyi(x) + 022/2(2;) + / U{x, t)u{t)dt, 
and choosing 

2qx^ 2grV2 

2/1 (x) = y/xJ_i/3( ), y2(x) = VxJl/3[ ), 

we find 

yi(%2W-y'i(%2(t) = ^. 

Applying (PT|) we obtain 



... ^ . /^^. .2qt3/2 9g:j3/2 ,2ga;3/\ ,2gt3/2 \ 

f/(x,t) - 2n^- (^J_i/3( JL_)Ji/3(Ji-_) - J_^/3(Jl__)J^/3(JL_)j 

Hence 



U{x,t)u{t)dt 




<q-\-^'^ / \u{t)\t-^'Ut, 
Jo 



and the result follows. D 

To find the constant of integration we use the value of J^{x) at two points: x — i^ 
and at the first zero x = ji,i. The following bounds are known [16] (see [3] for a 
review of recent results), 

(30) j-^^ = ^ + 2-1/3^,^1/3 + 02 3:2_2^^_,/3^ ^^^^ 

where Ug is s*'' positive zero of the Airy function Ai{—x). 
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In the sequel we set 7 = 2^^l'^a\ — 1.855757... Thus, we can write 

(31) U = v + iv^i^^e^^-^v-^i^. 

Lemma 12. 

(32) j^(^ + V/3)< 7^ ^>1 

\Sv 2 

Proof. Let j — j„i, by |3T|) we have 

J,(. + 7.V3) ^ _02 ^ ,-1/3 J^(, + ^,1/3 + ^2 ^ ,-1/3), 

Settmg 



by (fTS]) and 7 > -^i/F" we obtain 






where 

</)(e, S) = ^^ (4(8.5^ - 3) + 16{3S^ - 2)e + 4(6(53 _ 7)e^ + A{d'^ - 'S)e^ - 3e^)"'^^ 

Notice that (f){e, S) < 0(e, 7) since ^ > 7 and 

^0(.,^) = ^^0^(e,5)<O. 



Moreover, 



6„-y^ = 71.-2/3 + ^ ,-4/3 < 6^ 



whereas 

^</)(e,7)<0, 0<e<7. 
Thus, (/)(e,5) < (f>{0,j), and we get 

and the result follows. D 

Applying Lemma [11] to Ju we obtain 
Theorem 13. Let v > ^ , and let x = v + v^/^z, then 

,33) ..M = ^|^A.(-2./3.)..f°;"";f''. 

Vj/2/3 + 2 2i/2/3Vj/2/3 + 2 

Proof. Consider the function 

y{z) = \/v + i/i/3z J^iv + i^^/^z), 
which, as easy to check, satisfies the differential equation 

„^ ^ o . ^ 8z3 + 12iy2/322 _ 1 
y (z) + 2zy(z) = 4(^2/3 + ^)2 ^ W' 
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Lemmas [IT] and [5] yield 

(34) j/(z) = Vi(ciJ_i/3(V2C)+C2Ji/3(V2C)) +0i?(z), 

where C = ^ — , and 



z'^'R{z)= ^ 



4.V2J0 



8t^ + 12j/2/3t2 - 1 

(i.2/3 + t)2tl/4 y^^> 



dt < 



T7^-(^ -^) ' y — dt 



4V2 Jo (y'^'^ + i)3/2ii/4 

Z/V6 /- |8i3 + 12^.2/3^2 „ 1 



2v^ Jo iy'^''^ + 03/2(41.2/3^2 4^ 8^4/3i + 1)1/4^1/ 



dt < 



^ r r (3t/2/3 + 2t)t3/2dt /•" d±_ 

f^'^'^^Wo {iy^/^ + t)^/H2l,V3+ty/i^J^ 4Vt (1^2/3 +t)3/2 (21.2/3 +i) 1/4 /"< 

1 r r (3i/2/3 + 2t)i3/2dt [^ dt 

max ■ 



V^ [Jo (i.2/3+f)3/2(2j.2/3 + i)l/4'_yg 4-21/4^,7/6^ 

For i > the function 

3t/2/3 + 2t 

(1/2/3 +i)3/2(2i.2/3_^i)l/4 

is decreasing in t and therefore is less than 3 • 2^^^'^v^^/'^ , its value at i = 0. Thus, 
we obtain 

1 ( f^ 3i'^/2 /z 

(35) i?(z) < ^=-^ max jj^ ^^^^ dt, 



1/4. a; ' 25/4j,7/6 



max<^ - z^^'^,1^ ^/3 



25/4 0Fi; ]^ 5 

It is left to find the constants ci,C2. For x — u, that is for z = 0, we have 
R{z) — 0, and comparing ([M)) with PH)) we get 

2/(0) = Inn Vl(ciJ_i/3(V2C) +C2Ji/3(y2C)) = ^i 1™, V^ ' ^1/6^2 )^ = 
3i/3ci ^^,, 21/307 



2l/6r(|) ^ "^ ' 32/3r(|)(l/ + 02a)l/3' 

hence 

^ V2^ ^ V2v^/^ / _ 2^ 

^^ 3(z/ + 6i2a)i/3 3 V 3i/. 

Thus, using this and setting C3 = C2 — Ci, we can write 

(36) y{z) = 21/3^.1/6 v4z(-2i/3z)- 

^^ Ai{-2''^z) + C3 Vi Ji/3(V2 C) + e2i?(z) := 2'/^v''^ Ai{-2''\) + 7^. 
Now Ai(— 21/^7) — and by Lemma [T2l we have 



w ' 3 20 3vi^ 

yielding 

10 ,^ /- ^7T, 40 

C3 < -^ (2 + VI + iv-^/^) < 



3v^ 3a/j7 
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On the other hand, 

giving C3 > -7/v^. 
Inequahty (^01 ) gives 



7-27/12(4^ + 301/3)1/4 ' 

and using ^ and p^ we have 

VI|Ji/3(V2C)| < -M^^^ , 2V4vlzV4} • 

Combining these estimates with psp after some straightforward calculations one 
finds 

23max{M^ 
' ' "^ 2v^ 

and 1331 follows. D 

6. Sharper asymptotics 

The classical asymptotic given by (l?fl) does not makes much sense for x — 0{fj.) 
when the main term and the error are of the same order. Here using formula ([7]) 
we derive a different asymptotic expression with much smaller error term. It also 
leads to very sharp approximation of the Airy function Ai{—x). 

We'll need a few lemmas given in 18^. 

Lemma 14. Let f{x) satisfy the differential equation 

r'ix) + b\x)f{x)^o, 

where b{x) > and b"{x) exists on an interval T. 

Let g{x) = y/b(x)f{x), 

then for x £ I, provided the integral exists, 

3b'^{t)-2b{t)b"{t) 
463 (t) 

where B{x) = J b(t)dt, a £ I is arbitrary and \9\ < 1. 

Proof. Observe that g{x) satisfies the equation 

(38) g"-jg' + g{b'-e)^0, e ^ e{x) ^ — . 

The solution of the corresponding homogeneous equation 

is 

go — ci sinS(x) + C2 cosB{x). 

Solving formally p8p as a nonhomogeneous equation with the right hand side 
e{x)g(x) we get 

g{x) = go(x) + j ''-^^ ^ e{t)g{t)dt = 



(37) .g(a;) = ci cosS(x) + C2 sinS(a;) 



9{x) 



dt. 
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go{x) + 
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sm{B{x)-B{t)) 



m 



it)9{t) 



3&'2(i) - 2b{t)b"{t) 



4fe3(i) 



'W 



dt = 
dt. 



The normal form of differential equation ^ is 

/" + (1 - ^^-^)/ = 0, f^yGMx). 



Thus, for X > 1 /max{0, v"^ — j} we have 



6(a;) = 
\/x^+7i + y/]^ In 



2-2 _ j^2 _|_ 1. 



B(x) = 



^+^/|J.+x^ 



and 



y/x^~~jL + y/Jl arcsin ^ , 



1^1 < 



X > 



Theorem 15. For |zy| < i anrf a; > 0, 



(39) J^{x) ^J-ix'+ ti)~''^ cos {B{x) - uj^) + 

i^or \v\ > i and x > fi, 



12-KX (x2 +/i)3/2 



13/i 



(40) J,(x) = J-(a;2-/,)-V4eos(S(a;)-c.,)+L 

V TT 12V27r (a;2 - iifl'^ 

Proof. Since | Ji.(a;)| < J^ for |z^| < i , by dS?]) we have 

6z^ + ^ 



g(a;) = .90 + ^ 



M 



78^7, z3/2 (z2 + /i)9/4 

Comparing this with the standard asymptotic 

TT 

for large x one finds 

'2 



dz = go 



(41) 



/(x) = ^/xJ^{x) 



cos(a;-a;i.) +0(x"^/^), 



Cl 



smw^, C2 



COSWj^, 



yielding (jM)) . 

Similarly, for j/ > i and x > /i, using p7|) instead of (fT3|) . we obtain 



D 



9{x) =.90- 



8^ 7, (z2 - ^)5/2 



6z — /i 1 

^7:r dz = .90 + ^ 



Stt 



3x2 ^_ 2^ arcsin -^ 



3(x2 - ^)3/2 ^ 
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It is easy to check that the function 



3a;2 + 2^ arcsin ^ 



3(2;2 - ^)3/2 ^ 

is decreasing and positive. Therefore, by arcsin— > — , x > 0, we get 
3x2 ^ 2fi arcsin ^ 3x^ + 2^ 1 _ 



3(x2 - /i)3/2 ^ - 3(a;2 _ ^)3/2 a; 



13Ai (x + 2^x2 - A* )( Va;^ - fi-xf 13Ai 



6(x2 - /i)3/2 2x(x2 - ^)3/2 6(x2 - /i)3/2 ' 

and 

0(X) = On + t^ z= — . 

^^ ' ^ 12^/2^ (x2-/i)3/2 

Comparing this with the asymptotic for large x one finds 

ci — \ — smwy, Ci = -i/— coswy, 

V TT V TT 

and dlOl) foUows. D 



Remark 2. T/ie approximation given by ^JW *s strong enough to be matched with 
the solution in the transition region, as for x — v + v^'^ z the error term is of 
order 0{v~^'^ z~'''^). This allows one to find the constants of integrations ci, C2 in 
Theorem \13\ avoiding delicate claims like Lemma \JM We omit the details. 

Similarly to Corollary [T] we obtain 

Corollary 2. 

(42) Ai{-x) = 



v^ (16x3 + 5)1/4 0Fxi/4(16x3 + 5)3/2 ' 

provided x > 0. 

The argument of the cosine in (|5^ and (|42l) (but not in (HH])) can be simplified 
at the cost of a weaker numerical constant at the error term. Namely, it is easy to 
verify the following elementary inequalities: 

2 

y/x^ + ^i + ^ In ^ "", J =^"^ + ^'^' ^>0, 

y/JJ+y/fi + X^ 2x 24x3 



Vl6x3 + 5 VS , Vl6x3 + 5 + V5 _2 3/2 _ _5_ -3/2 n2^^ -9/2 
6 6 4x3/2 3 48 9216 

Since | cos(x + e) — cosx| < e, we obtain 

^ ' ' Vtt (x2+^)i/4 24V2^x3(x2 + /i)i/4' '11-2 

2^/icos(-x3/2 - A2:-3/2 _ IL\ K 
(44) Ai{~x) = ^ '^3 4^^ 1^ ' « ^ 



0F(16x3 + 5)1/4 9V7Fx4(16x3 + 5)1/4 
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In particular, the last formula yields rather sharp approximations for the zeros 
of Airy (e.g. already for the first zero the error is less than 0.00122), and, in view 
of the inequality dSH]) . the Bessel function. 

Theorem 16. 

where as is s*'' positive zero of Ai{—x) and m = (12s — 3)7r. 

Proof. Elementary arguments show that |sina;| < e implies sin(x + ^^) — for 
some 6, \9\ < 1. Therefore, Ai{—x) = means 

I ^'^' - Js ^"'^' - T + ^ S ^"'^' = "'' ' = °' ^' ^' - 

Since ai — 2.33..., we may assume x > 2. Then 

36 



Aa;-3/2 > ^3,-9/2 



hence 

Z 0/9 OTT 

_a;J/- > _+7rs, 
and we obtain the following estimate 

^^-9/2^ 640 



36 2437r2(4s + 3)3 ■ 

Thus, for (s + 1)*'* zero this gives the equation 

2^/9 5 o/n Stt ^ 6407r 37r 

- x''/'^ x"-^'-^ = h TTs + 61 :^ h ITS + e, 

3 48 4 9m3 4 ' 

with the relevant solution 

X = 16"^/^ fm + 12e + ^(m + 12e)2 +40] 
After some calculations one gets 

a; = 16"^/^ (m + Vm^ + 4o) (l + e^=^^] = 
V / V Vm2 + 40 / 

>2/3 „ 12807r 



16"^/^ [771+ ^^2+40] +( 



9TO3(m2+ 40)1/6 ■ 
This completes the proof. D 



Formula p5|) can be simplified at the cost of slightly weaker numerical constant. 
Namely, as one can checks 



< -(^2 + 20)^/'"^ - 16^2/^ (m + Vm^ + 40 ] 



2/3 25 

< 



3to3(to2+ 40)1/6 
yielding 



4^ m3(m2 + 40)1/6 ' 

Finally, comparing the numerical values of the zeros of Ai(—x) with (P5|) leads to 
the following conjecture: 

as < 16-2/3 (m + v/m2 + 4o] . 
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